Topological code Autotune 



Austin G. Fowler 1 , Adam C. Whiteside 1 , Angus L. Mclnnes 1 , Alimohammad Rabbani 2 
1 Centre for Quantum Computation and Communication Technology, 
School of Physics, The University of Melbourne, Victoria 3010, Australia, 
2 Sharif University of Technology, Department of Computer Engineering, Tehran, Iran 

(Dated: February 29, 2012) 

Tailoring a fault-tolerant quantum error correction scheme to a specific physical architecture can 
be a laborious task. We describe a tool Autotune capable of analyzing and optimizing the classical 
processing for an arbitrary 2-D qubit architecture making use of arbitrary circuits implementing 
either the surface code or progressively generated slices of a 3-D topological cluster state with 
arbitrary stochastic error models for each quantum gate. Autotune is designed to facilitate precise 
study of the performance of real hardware running topological quantum error correction. 



Large-scale fault-tolerant simulations of the surface 
code indicate gate error rates between 0.2% and 0.5% are 
sufficiently low to enable practical overhead, high reliabil- 
ity quantum computation [T]. This is tantalisingly close 
to experimentally achieved two-qubit gate error rates of 
2% [2], the best achieved to date in a system with the 
potential to implement the required 2-D array of qubits. 
This motivates the serious study of mapping topologi- 
cal quantum error correction (TQEC) schemes to physi- 
cal hardware to enable realistic engineering trade-offs to 
be determined and optimizations found. We would very 
much like to collaborate with any experimentalist with a 
potentially 2-D qubit system and an interest in making 
use of TQEC. Given the transversely invariant nature of 
TQEC, experiments involving as few as two qubits can 
be sufficient to determine whether a system could suc- 
cessfully implement TQEC. 

This paper is not self-contained. We assume the reader 
has basic familiarity with surface code quantum compu- 
tation and topological cluster states [3-6 . Our primary 
focus will be explaining how our tool Autotune works and 
is used. 

The discussion is organized as follows. In Section |TJ 
we describe how error models are defined for each quan- 
tum gate and the current range of phenomena they can 
encompass. In Section [TTl we describe how errors are 



tracked during simulation. In Section III the concept of 
a set of measurements is introduced, with a -1 product 
of measurements in a set indicating a detection event, 
defined in Section [IV| Away from boundaries, every sin- 
gle error leads to a pair of detection events. We store 
this information in detail in a nest (Section [V]) . This de- 
tailed information can be distilled on the fly into a lattice 
(Section VI), however this is somewhat computationally 
expensive (a large constant overhead) and so we describe 
how to distil this information once during a bootup phase 
into a recipe (Section VII[ ) that enables high-speed con- 
struction of the appropriate lattice for matching without 
the expensive tracking of errors. Having described Au- 
totune, we give examples of its use to analyze the stan- 
dard surface code (Section Villi and progressively gen- 
crated 3-D topological cluster states (Section IX). Sec- 
tion [X] concludes with a discussion of our planned future 
work. 



I. ERROR MODELS 

Autotune is capable of handling any error model with 
outcomes that can be described by a single integer e per 
qubit. For example, I = 0, X = 1, Z = 2, Y = 3, 
leaked to a non-computational state = 4, qubit lost = 
5. Only single-qubit and two-qubit gate error models 
are currently supported, however this could easily be ex- 
tended. The user can specify how e is transformed by 
each gate, for example controlling whether a CNOT be- 
tween a leaked qubit and a non-leaked qubit results in 
two leaked qubits, or no effect on the non-leaked qubit, 
or any other effect describable by single integers. It must 
be possible to deterministically say whether each error e 
contains X and/or Z or neither. Y contains both, I nei- 
ther, and a leaked qubit may look like either or neither 
depending on the underlying physics. 

An example of a Pauli channel CNOT error model file 
is shown below. 
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The first line states the number of qubits n q the gate is 
applied to. The second line states the value x the relative 
strengths Sj of the various errors should be normalized to 
sum to. This makes it easy to handle different gates with 
different overall probabilities of error, e.g. x\ q = X2 q /10. 
The third line states the number of different errors n e 
in the model. The next n e lines contain n q + 1 integers 
specifying the relative strength of that error and the value 
of e to apply to each qubit. The user can specify exactly 
how different errors combine e\e2 — e^. Each time a 
quantum gate is called, it is passed an error rate p and an 
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error of any kind is applied with probability px. Error i is 
then applied with mutually exclusive probability Sj/Ss,. 
The final line is the gate duration in arbitrary units. 

II. ERROR TRACKING 

Error models are not only used to generate stochastic 
errors, we also deterministically track a representative set 
of errors. Specifically, given an arbitrary single-qubit er- 
ror model we sum the relative probabilities of errors that 
contain X and those that contain Z and track these two 
non mutually exclusive errors only. Given an arbitrary 
two-qubit error model we track the non mutually exclu- 
sive errors IX, XI, XX, IZ, ZI and ZZ. This is done 
simply to keep the number of tracked errors under con- 
trol. Transforming the CNOT error model in Section [T] 
in this manner yields 
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Every time a gate is applied, all nonzero relative proba- 
bility pure X and pure Z errors are generated and added 
to the list of errors on each qubit touched by the gate. 
Relative probabilities are stored in each error scaled by 
the p passed to the gate. Each error is given a unique la- 
bel. Multiple-qubit errors are represented by single-qubit 
errors on each qubit each with the same label. If one ex- 
ecutes a long sequence of unitary gates, the number of 
errors per qubit that need to be tracked will grow without 
bound. Each unitary gate transforms all errors present 
on all of the qubits it touches. An S gate will transform 
X errors into Y errors and vice versa. Multiple-qubit 
gates can create new propagated errors which will have 
the same label as the original and can combine or cancel 
multiple errors on a single qubit with the same label. For 
example, CNOT((7i, q2, p = 0.01) applied to qubits 

qi -> (X, 0.005, 0) 

^ (Z, 0.00238,1) (3) 

q 2 -> (Y, 0.00238, 1) (4) 

where (A, p sr , L) represents the error type, scaled rel- 
ative probability and label and the arrows represent a 
linked list, will result in 

qi -> (X, 0.005, 0) 

^ {X, 0.005, 3) 

^ {Z, 0.00373, 5) 

^ {Z, 0.00238, 6) (5) 



q 2 -> (X, 0.005, 0) 

^ (Y, 0.00238, 1) 

^ (X, 0.00385, 2) 

^ (Z, 0.00332, 4) 

^ (Z, 0.00238, 6) (6) 



III. SETS OF MEASUREMENTS 

Autotunc currently supports only single-qubit mea- 
surements in the X and Z bases. Mx applied to eq. [6] 
will create a measurement 

m ->■ (Z, 0.00238, 1) 
(Z, 0.00332, 4) 

(Z, 0.00238, 6) (7) 

The X errors and X components of Y errors have been 
removed. All errors will be removed from the qubit. To 
use the qubit again it must be explicitly initialized. If the 
qubit does not need to be used immediately after being 
measured, Autotune provides a dead gate that advances 
the qubit in time but does not generate or track any 
errors. This models incoherent evolution. 

Every measurement is associated with either two sets 
or a single set and a boundary. A set of measurements 
has the property that the product of the measurement 
results (+1 or -1) indicates whether a chain of errors has 
ended nearby. In the standard surface code, sets contain 
consecutive pairs of syndrome qubit measurements. In a 
3-D topological cluster state, sets contain the measure- 
ments on the faces of individual primal and dual cells. 
See Fig. [T] Sets must be specified by the user. 

Sets can also be associated with boundaries. The bot- 
tom layer of sets in Fig. [T^, is associated with the primal 
initial time boundary. In the second layer of sets, sets in 
the left row are associated with the left primal boundary, 
those in the right row are associated with the right primal 
boundary. We use the terminology primal/dual instead 
of rough/smooth as used in [5] to ensure uniform termi- 
nology when discussing both the surface code and 3-D 
topological cluster states. Note that the middle row of 
sets in the second layer is not associated with any bound- 
ary. The association of sets with boundaries is currently 
manually user specified. 



IV. DETECTION EVENTS 

When all measurements in a set have been performed, 
further processing is triggered. A measurement may con- 
tain many errors. A set may contain many measure- 
ments. Autotune determines which errors with the same 
label appear an odd number of times. For each such la- 
bel a detection event is generated. Fig. [2^, contains an 
example of a single error leading to a pair of detection 
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FIG. 1: a.) 2-D surface code (grey). Time runs vertically. In 
the front left corner, the transversely invariant initialization, 
4 CNOT, then measurement pattern is shown. For clarity, 
only the measurement gate of this pattern is shown elsewhere 
in space-time. Long ovals encircle pairs of measurements and 
represent sets. Arrows indicate associations with boundaries. 
The bottom layer of sets is associated with the initial time 
boundary. The second layer left and right rows of sets are 
associated with the left and right spatial boundaries respec- 
tively, b.) A single 3-D topological cluster state cell. Sets 
contain six measurements away from boundaries of the lat- 
tice. 



events. Detection events are stored in a hash table to en- 
able one to quickly determine whether a detection event 
with a given label has already been generated. Pairs 
of detection events immediately trigger the creation of 
sticks, described in the following Section. 

Errors near a boundary can lead to single detection 
events (Fig. [2]d) . Such single detection events must gen- 




FIG. 2: 2-D surface code (grey). Time runs vertically. 
Squares represent initialization to |0), circles represent Z basis 
measurement, a.) A single error leading to a pair of detection 
events, b.) An error leading to a single detection event due 
to proximity to a boundary of the lattice. 

erate sticks leading to the nearby boundary. One must 
decide with care when to conclude that a detection event 
is unique and that no matching detection event will be 
generated in the future. To deal with this, wc define 
a measure of error detection progress big_t to increment 
only when every stabilizer of the code has been measured 
at least once. An example of an error leading to detec- 
tion events two big_t in the future is shown in Fig. [3] 
This is the maximum delay possible. Detection events 
three big_t in the past that are unique are guaranteed to 
remain so. 

Detection events can be primal or dual depending 
whether they arc associated with primal or dual sets 
which are in turn associated with primal or dual stabilizer 
measurements. It is a good idea to design stabilizer mea- 
surement circuits with the property that single errors do 
not lead to the creation of more than one pair of primal 
and one pair of dual detection events. This ensures that 
minimum weight perfect matching [7J [5] is well-suited to 
correcting errors generated during the execution of the 
circuits. 



V. NESTS OF BALLS AND STICKS 

When a set is processed to generate detection events, 
it is said to be finalized. At this point in time, a ball 
is generated and associated with the finalized set. This 
ball represents a location in space-time. A pair of detec- 
tion events leads to the creation of a stick between its 
associated balls. A stick represents a potential connec- 
tion between a pair of space-time locations arising from 
a single error. Many single errors can lead to the same 
stick [5]. The probability of a given stick, given a list of 
errors leading to it each with its own probability (ej,Pi), 
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FIG. 3: Example of a single error leading to two detection 
events two big_t in the future. A detection event has a big_t 
equal to that of its latest measurement. big_t increments after 
the every stabilizers has been measured at least once (circles 
represent measurement, squares represent initialisation). 



is therefore, to first order 



Pstick = "^2 Pi Yii 1 ~Pj)- 



(8) 



We have a simple Blender based visualization tool for 
nests, a screen shot of which can be found in Fig. [4] 
A nest contains full tracking information of which gate 
led to which collection of errors and which of those er- 
rors led to which sticks. This is highly useful, however 
computationally cumbersome. Autotune does delete all 
errors, measurements, sets, detection events, balls and 
sticks when they are no longer required, which keeps total 
memory required finite, however it remains challenging 
to generate the nest fast enough to obtain good statistics 
in simulations, let alone keep pace with a real quantum 
computer. 



VI. LATTICES OF DOTS AND LINES 

Motivated by the difficulty of rapidly generating nests, 
we first cut down the data stored therein to the mini- 
mum required by the minimum weight perfect matching 
algorithm. For every ball in a nest we create a dot, which 
again simply represents a space-time location. For every 
stick, we create a line. Lines connect the dots corre- 
sponding to the balls the stick connected, however they 
contain just one number, the weight w = — ln(p s tick)' 

The lattice is kept for much longer than the nest as 
with low probability the matching algorithm can require 



all prior matching history to correctly match the latest 
data. A sufficiently long history must be kept to ensure 
the probability of requiring more is negligible. This is 
possible as the probability of requiring additional data in 
the past decreases exponentially with kept history size. 
The simplicity of the lattice keeps the memory required 
low. 



VII. RECIPES OF OFFSETS, 
LAYERS 



BLOCKS AND 



The visualized nest in Fig. [4] has a great deal of reg- 
ularity in its structure. This can be exploited to enable 
direct generation of lattices, avoiding the need to gener- 
ate expensive nests on the fly. During a boot-up phase, 
Autotune analyzes each new stick and stores new unique 
sticks as an offset that contains only the geometric in- 
formation associated with the stick — just enough to 
create a line. Furthermore, each new ball is analyzed 
to determine if it corresponds to a new pattern of off- 
sets. Such unique patterns are stored as blocks of offsets. 
Finally, each round of error detection is analyzed with 
structurally unique rounds stored as layers of blocks. All 
data is then stored in a recipe that contains all necessary 
information to rapidly create any part of the lattice. 

Currently, Autotune is capable of analyzing either er- 
ror detection circuitry that eventually leads to identi- 
cal repeated layers or a finite number of rounds of er- 
ror detection that can have any structure whatsoever. 
The former could easily be extended to a finite number 
of cyclically repeated layers, however it remains unclear 
whether one could avoid generating full nests when simu- 
lating probabilistic error detection in which all syndrome 
measurements can take a randomly variable amount of 
time. 



VIII. EXAMPLE: THE SURFACE CODE 

Consider the stabilizer measurement circuits shown in 
Fig. [5] These circuits are highly physical, with initializa- 
tion to |0) and Z basis measurement only, and specifically 
no |+) initialization or X basis measurement as was used 
in [TJ [TU] . Note the use of incoherent evolution dead gates 
D to ensure that Z-stabilizer measurement takes as long 
as A"-stabilizer measurement. We do not assume that 
all gates take equal amounts of time. Gate D is triv- 
ially made to take the same time as H, however all other 
types of gates are allowed to take arbitrary amounts of 
time. This implies we need four different duration iden- 
tity gates — Tjnit, Iff, Imcas and I C not- These different 
identity gates are applied to the data qubits at appro- 
priate times while syndrome qubits are undergoing the 
corresponding operations. 

Assuming equal probability p depolarizing noise on all 
unitary gates and incorrect initialization and measure- 
ment with probability p, the low distance performance is 
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FIG. 4: 3-D Blender model of the probability of sticks. Thicker cylinders represent more probable sticks. Time runs vertically, 
the (a;, y) — (0, 0) position is the back corner. 
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FIG. 5: Physically reasonable syndrome measurement cir- 
cuits. D gates represent incoherent evolution. 



shown in Fig. [6] We expect a significantly lower threshold 
error rate than that observed in [JJ due to the additional 
gates present in the syndrome measurement circuits. The 
high distance performance about the threshold error rate 
is shown in Fig. [7] We observe a threshold error rate of 
approximately 0.57%. 



IX. EXAMPLE: TOPOLOGICAL CLUSTER 
STATES 



A 3-D topological cluster state can be progressively 
generated using two layers of qubits and a repeating 6- 
step cycle of measurements and interactions. This is de- 
scribed in detail in Fig. [8] A single layer of qubits is 
sufficient provided one can have some overlapping but 
still nearest neighbor interactions (Fig. |9| . Sets require 
more care to construct as they now contain up to six mea- 
surements, one on each face of each primal and dual cell. 
The generic form of the primal and dual cell structure of 
the cluster state we will simulate is shown in Fig. [TUJ 

We assume the same error model as that used for the 
surface code in the previous Section, however we per- 
mit the use of initialization to |+) and JWj. Initializa- 
tion to |+) produces |— ) with probability p. Mx reports 
the wrong cigenstate with probability p. Measurement 
must be followed by initialization. The dual nest result- 
ing from the described gate sequence and error model for 
a distance 4 topological cluster state is shown in Fig. 11 
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FIG. 6: Probability of logical X error per round of error correction for different square surface code distances d. Quadratic 
and cubic asymptotic curves are shown for distances 3 and 5 respectively. 
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FIG. 7: Probability of logical X error per round of error correction for large distances around the threshold error rate of 
approximately 0.57% 



Note that only one type of logical error is now possible as 
we have a single unbroken U-shaped primal boundary for 
all time and the only nontrivial errors connect the two 
opposing dual boundaries. 

As was done for the surface code, for low error rates 
we can compare the actual performance to the expected 
asymptotic behavior. A logical error in a distance 3 code 
requires a minimum of two physical errors to occur, so 



for low error rates a distance 3 code is expected to have 
an error rate quadratic in the depolarizing probability 
p. Similarly, distance 5 codes should converge to a cu- 
bic curve and distance 7 to a quartic curve. It can be 



seen from Fig. 12 that our results converge well to the 
expected single-term polynomials. 

In addition to the asymptotic behavior, the thresh- 



old error rate is of interest. Fig. 13 shows the behavior 



FIG. 8: A unit cell of the interaction pattern of a double 
layer of qubits generating a 3-D cluster state. Small black 
dots correspond to qubits. Most of the qubits in the second 
layer have been omitted for clarity. The numbers at the top 
right and bottom left of each qubit indicate the time steps in 
which it is initialized and measured, respectively. Each line 
corresponds to a Cz gate and is labeled with the time step(s) 
in which this gate is applied. The other layer comprises an 
identical ordering of qubits and gates, but shifted in time by 
three steps (such that a gate labeled 2 is executed in time 
step 5 and vice versa, etc.) and shifted north, east and up as 
indicated by the arrow. Upward Cz gates in the first layer 
become downward after shifting. 
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FIG. 9: A single-layer implementation of the two layers used 
to progressively generating a 3-D topological cluster state. 
Note that while some interactions cross, they are still techni- 
cally nearest neighbor. 



around threshold for the error model described above. 
The observed threshold area of approximately 0.58% is 
lower than the threshold obtained in [3] due to the more 
severe error model. Fig. [14] is a repeat of the same cal- 
culation using depolarizing noise after initialization and 
before measurement leading to an approximate threshold 
error rate of 0.66%, in good agreement with the 0.67% 
value in eq. 23 of 0]. 




FIG. 10: Primal cells (top) and dual cells (bottom) in a 
distance 4 example of the topological 3-D cluster states we 
simulate. Dual boundaries are present on the left and right 
sides and a single U-shaped primal boundary covers the bot- 
tom, front and back. Time runs vertically. The top temporal 
surface should not be considered a boundary as the cluster 
state is being extended in this direction. 



X. CONCLUSION 

We have described a code library Autotune that is ca- 
pable of handling in a natural manner both fully fault- 
tolerant surface codes and 3-D topological cluster states. 
Generality is achieved through the definition of sets of 
measurements which are used to detect error chain end- 
points. Arbitrary syndrome measurement circuits are 
supported along with arbitrary stochastic error models 
for each gate. The code is highly efficient, with runtime 
comparable to that reported in 

In future work, we plan to take correlations between 
errors into account, parallelize the core matching en- 
gine, and develop the capability to analyze complex fault- 
tolerant circuits consisting of many braided defects. 
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FIG. 11: The dual nest associated with the cluster state of Fig. [l0| 
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FIG. 12: Logical error rate per layer of cluster state cells for distances 3, 5, 7 and 9. Dashed lines are polynomials showing 
the expected quadratic, cubic and quartic single-term polynomial asymptotic behaviors. 
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FIG. 13: Logical error rate per layer of cluster state cells for a large range of distances around the threshold error rate. 
Initialization and measurement fail with probability p. The threshold error rate for this error model is approximately 0.58%. 
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FIG. 14: Logical error rate per layer of cluster state cells for a large range of distances around the threshold error rate. 
Initialization and measurement suffer depolarizing noise with probability p. The threshold error rate for this error model is 
approximately 0.66%. 
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